Search for: All records

Abstract This article studies two particular algorithms, a relaxation least squares algorithm and a relaxation Newton iteration scheme, for reconstructing unknown parameters in dissipative dynamical systems. Both algorithms are based on a continuous data assimilation (CDA) algorithm for state reconstruction of Azouaniet al(2014J. Nonlinear Sci.24277–304). Due to the CDA origins of these parameter recovery algorithms, these schemes provide on-the-fly reconstruction, that is, as data is collected, of unknown state and parameters simultaneously. It is shown how both algorithms give way to a robust general framework for simultaneous state and parameter estimation. In particular, we develop a general theory, applicable to a large class of dissipative dynamical systems, which identifies structural and algorithmic conditions under which the proposed algorithms achieve reconstruction of the true parameters. The algorithms are implemented on a high-dimensional two-layer Lorenz 96 model, where the theoretical conditions of the general framework are explicitly verifiable. They are also implemented on the two-dimensional Rayleigh–Bénard convection system to demonstrate the applicability of the algorithms beyond the finite-dimensional setting. In each case, systematic numerical experiments are carried out probing the efficacy of the proposed algorithms, in addition to the apparent benefits and drawbacks between them. 

We study a class of semi-linear differential Volterra equations with polynomial-type potentials that incorporates the effects of memory while being subjected to random perturbations via an additive Gaussian noise. We show that for a broad class of non-linear potentials, the system always admits invariant probability measures. However, the presence of memory effects precludes access to compactness in a typical fashion. In this paper, this obstacle is overcome by introducing functional spaces adapted to the memory kernels, thereby allowing one to recover compactness. Under the assumption of sufficiently smooth noise, it is then shown that the statistically stationary states possess higher-order regularity properties dictated by the structure of the nonlinearity. This is established through a control argument that asymptotically transfers regularity onto the solution by exploiting the underlying Lyapunov structure of the system in a novel way. 

In this article, some elementary observations are made regarding the behavior of solutions to the two-dimensional curl-free Burgers equation which suggest the distinguished role played by the scalar divergence field in determining the dynamics of the solution. These observations inspire a new divergence-based regularity con- dition for the two-dimensional Kuramoto–Sivashinsky equation (KSE) that provides conceptual clarity to the nature of the potential blow-up mechanism for this system. The relation of this regularity criterion to the Ladyzhenskaya–Prodi–Serrin-type cri- terion for the KSE is also established, thus providing the basis for the development of an alternative framework of regularity criterion for this equation based solely on the low-mode behavior of its solutions. The article concludes by applying these ideas to identify a conceptually simple modification of KSE that yields globally regular solu- tions, as well as providing a straightforward verification of this regularity criterion to establish global regularity of solutions to the 2D Burgers–Sivashinsky equation. The proofs are direct, elementary, and concise. 

This article is concerned with the problem of determining an unknown source of non-potential, external time-dependent perturbations of an incompressible fluid from large-scale observations on the flow field. A relaxation-based approach is proposed for accomplishing this, which makes use of a nonlinear property of the equations of motions to asymptotically enslave small scales to large scales. In particular, an algorithm is introduced that systematically produces approximations of the flow field on the unobserved scales in order to generate an approximation to the unknown force; the process is then repeated to generate an improved approximation of the unobserved scales, and so on. A mathematical proof of convergence of this algorithm is established in the context of the two-dimensional Navier–Stokes equations with periodic boundary conditions under the assumption that the force belongs to the observational subspace of phase space; at each stage in the algorithm, it is shown that the model error, represented as the difference between the approximating and true force, asymptotically decreases to zero in a geometric fashion provided that sufficiently many scales are observed and certain parameters of the algorithm are appropriately tuned. 

Motivated by recent progress in data assimilation, we develop an algorithm to dynamically learn the parameters of a chaotic system from partial observations. Under reasonable assumptions, we supply a rigorous analytical proof that guarantees the convergence of this algorithm to the true parameter values when the system in question is the classic three-dimensional Lorenz system. Such a result appears to be the first of its kind for dynamical parameter estimation of nonlinear systems. Computationally, we demonstrate the efficacy of this algorithm on the Lorenz system by recovering any proper subset of the three non-dimensional parameters of the system, so long as a corresponding subset of the state is observable. We moreover probe the limitations of the algorithm by identifying dynamical regimes under which certain parameters cannot be effectively inferred having only observed certain state variables. In such cases, modifications to the algorithm are proposed that ultimately result in recovery of the parameter. Lastly, computational evidence is provided that supports the efficacy of the algorithm well beyond the hypotheses specified by the theorem, including in the presence of noisy observations, stochastic forcing, and the case where the observations are discrete and sparse in time. 

This paper studies a family of generalized surface quasi-geostrophic (SQG) equations for an active scalar \begin{document}$$ \theta $$\end{document} on the whole plane whose velocities have been mildly regularized, for instance, logarithmically. The well-posedness of these regularized models in borderline Sobolev regularity have previously been studied by D. Chae and J. Wu when the velocity \begin{document}$ u $$\end{document} is of lower singularity, i.e., \begin{document}$$ u = -\nabla^{\perp} \Lambda^{ \beta-2}p( \Lambda) \theta $$\end{document}, where \begin{document}$$ p $$\end{document} is a logarithmic smoothing operator and \begin{document}$$ \beta \in [0, 1] $$\end{document}. We complete this study by considering the more singular regime \begin{document}$$ \beta\in(1, 2) $$\end{document}$. The main tool is the identification of a suitable linearized system that preserves the underlying commutator structure for the original equation. We observe that this structure is ultimately crucial for obtaining continuity of the flow map. In particular, straightforward applications of previous methods for active transport equations fail to capture the more nuanced commutator structure of the equation in this more singular regime. The proposed linearized system nontrivially modifies the flux of the original system in such a way that it coincides with the original flux when evaluated along solutions of the original system. The requisite estimates are developed for this modified linear system to ensure its well-posedness.